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ABSTRACT 

I show how the existing observational data on Local Group dwarf galaxies can be 



CN ' used to estimate the average star formation law during the first 3 Gyr of the history 

. of the universe. I find that the observational data are consistent with the orthodox 

, Schmidt law with a star formation efficiency of about 4 percent if the star formation 

I is continuous (during the first 3 Gyr). The efficiency is proportionally higher if most 

of the gas in the dwarfs was consumed (and never replenished) in a short time interval 
well before the universe turned 3 Gyr. 

> ' 

Subject headings: cosmology: theory - galaxies: formation - stars: formation 
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O ■ 1. Introduction 

C^l Direct observational measurements of the star formation rate in the very first objects in 

O ' the universe is exceptionally difficult, and even with the recent observational and technological 

' advances, is well beyond reach. Nevertheless, it is possible to reconstruct the general features 

■ of the star formation at high redshift in much the same way as paleontologists reconstruct the 



prehistoric life from fossils. 



In the currently widely accepted hierarchical clustering paradigm the characteristic mass 
5^ I scale of cosmological objects increases with time: low mass objects form first, and more massive 

objects form later. Thus, a considerable fraction of the dwarf galaxies may be very old, and may 
represent the fossils of the very first star formation in the universe. Although these dwarf galaxies 
are small and faint, and thus hard to observe, a vast amount of detailed information exists on the 
Local Group dwarf galaxies, as summarized in several recent reviews (Grebel 1998; Mateo 1998; 
van den Bergh 1999). The first two reviews are particularly useful for my purpose, as they contain 
the star formation histories of a large fraction of all Local Group dwarfs. A careful investigation 
of the available data reveals a striking feature common to the vast majority (more than 90%) of 
all dwarf galaxies with known star formation histories: a sharp decline in the star formation rate 
about 10 Gyr ago. 



There are at least five possible explanations for such a feature in the star formation history. 
First, it may be attributed to the tidal effects from the Galaxy or Andromeda. However, in 
the hierarchically clustering universe neither the Galaxy nor Andromeda existed at that time 
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as entities, but rather as an interacting collection of proto-galactic clumps (see, for example, 
Contardo, Stcinmctz, & Fritze-von Alvensleben 1998), and thus they were not able to have a large 
impact on the star formation histories of almost all dwarf galaxies. 

A second possible explanation is the effect of supernovae, which could have expelled the gas 
from the low potential wells and quenched the star formation rate. However, it is hard to believe 
that the supernova activity was so well synchronized in 90% of all dwarf galaxies as to expel the 
gas almost simultaneously (this is elaborated further in the Conclusions). 

Since the drops in the star formation rates in dwarf galaxies are essentially simultaneous, 
it is plausible that this was triggered by a process that affected the whole universe. One such 
process is reionization. It is well established that photoheating of the intergalactic medium during 
cosmological reionization sharply reduces the gas fractions in the low mass objects (Thoul & 
Weinberg 1996; Quinn, Katz, & Efstathiou 1996; Weinberg, Hernquist, & Katz 1997; Navarro &i 
Steinmetz 1997; Gnedin 2000). This, in turn, will lead to the drop in the star formation rates in 
dwarf galaxies, including the Local Group dwarfs. This proposal however encounters a difficulty 
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because an age of 10 Gyr corresponds to about z = 2 (in a 13 Gyr old universe with h%' = 0.4), 
which is too late for cosmological reionization. However, if the drop in the star formation rate in 
dwarf galaxies occurred some twelve billion years ago, the reionization explanation could work. 
This requires that the measured stellar ages are off by about 2 Gyr, but the question of whether 
this is plausible or not is outside the scope of this paper. 

The fourth explanation is based on the observational evidence (which is admittedly not 
completely compelling) that the universe experienced a second epoch of reheating at about z = 3 
(Ricotti, Gnedin, & Shull 2000; Schaye et al. 2000). This second reheating (plausibly attributed 
to the reionization of helium) has a similar effect on dwarf galaxies as hydrogen reionization at 
z ~ 7 — 10, except that it takes place somewhere between z = 2 and z = 3. 

It is also possible that photoionization associated with the burst of star formation that made 
the bulges of the Galaxy and Andromeda ( "local reionization" ) inhibited the star formation in the 
surrounding dwarfs (van den Bergh 1994). While possible in principle, it is not clear whether this 
explanation is plausible within the hierarchical clustering paradigm, or whether it agrees with the 
observed metallicity distribution of the Galactic bulge (McWilliam & Rich 1994). 

Whatever the explanation of the observed simultaneous drop in the star formation rate in 
90% of Local Group dwarfs is, it is not essential to the measurement described in this paper. The 
only thing that matters is that there exists a subset of the Local Group dwarfs that had little 
star formation in the last 10 Gyr, so that their stellar content provides a fossil record of the early 
universe. 
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Table 1. Local Group Dwarfs Used in this Analysis 



Galaxy 




1 / T / / T '-{wYt 

\og{Lv/{Icri)T 


log(/cPc7L0)''''= 


t d 
/>10 


NGC 3109 


0.2-0.5'= 


1.55 ±0.23 


-1.75 ±0.20(0.33) 


0.7-0.5 " 


NGC 205 


0.2 


1.69 ±0.27 


-0.37 ±0.20(0.25) 


0.8 


NGC 147 


0.3 


1.84 ±0.15 


-0.41 ±0.12(0.17) 


0.7 


Ursa Minor 


0-0.2^ 


0.78 ±0.26 


-2.22 ±0.21(0.23) 


0.9 


Sextans 


0.1 


0.82 ±0.24 


-2.70 ±0.21(0.28) 


0.8-0.7 


Antlia 


0-0.2° 


0.95 ±0.12 


-1.80 ±0.10(0.12) 


0.9 


Sculptor 


0.1 


1.47 ±0.20 


-1.26 ±0.17(0.19) 


0.9 


And I 


0-o.r 


1.00 ±0.12 


-2.05 ±0.10(0.11) 


0.9 


And II 


0.2 


1.21 ±0.12 


-1.77 ±0.23(0.25) 


0.9 


And III 


0-0. r 


1.03 ±0.14 


-1.75 ±0.10(0.11) 


0.9 



'^The ratio of the star formation rate just after the drop at 10 ± 1 Gyr ago to the star formation 
rate just before the drop. Estimated from Fig. 8 of Mateo (1998) or Fig. 5-6 of Grebel (1998). This 
column is shown for illustration only, and is not used in the analysis. 

•^The data are taken from Tables 2-4 of Mateo (1998). 

■^The errors in parenthesis are corrected by a factor of and outside of parenthesis are uncorrected. 

"^The fraction of all stars formed more than 10 Gyr ago. Estimated from Fig. 8 of Mateo (1998) 
or Fig. 5-6 of Grebel (1998). 

''The first number is from Mateo (1998) and the second one is from Grebel (1998); for />io the 
arithmetic mean is used in the analysis. 

2. Results 

Because the star formation histories of the Local Group dwarfs are quite diverse, only a 
subset of all galaxies can be used, namely those that have a large fraction of their stars formed 
more than 10 Gyr ago. This limits the available sample to seventeen out of thirty two shown in 
Fig. 8 of Mateo (1998) and Fig. 5-6 of Grebel (1998). Out of a total of seventeen, seven (Fornax, 
Sextans B, NGC 185, Phoenix, GR 8, Draco, and Tucana) either have some of the essential data 
missing or Mateo (1998), Grebel (1998) and Hernandez, Gilmore, &: Valls-Gabaud (2000) strongly 
disagree on their star formation histories, so only the ten listed in Table |I| are used in the present 
analysis. The other seven can be added to this analysis when the missing or disgreed upon data 
become available. 



According to the Schmidt law, on sufficiently large scales the star formation rate can be 
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parametrized as a power-law of the gas density, 

^ - (1) 

where is the mass density of stars, and pg is the mass density of gas. Integrating equation 
(|l]) over the total volume of a galaxy and time, and assuming that the mass-to-light ratio is 
independent of radius, I can obtain the following equation connecting the total luminosity of a 
galaxy at some moment tf with the total baryon core density p^ ,, a-t the same time. 



^0 \Pb,c J \rcJ rc 

where is the luminosity of a dwarf galaxy, Tc is the core radius, and Atcfr is the effective time 
interval, Atgfr = Jo^ J PgdV dt/ J p'^{tf)dV, which incorporates all the details of the increase (due 
to infall) and decrease (due to consumption by star formation) of the gas. 

It is important that I use the core baryon density rather than core gas density. The gas in 
the center can be almost entirely consumed, whereas the baryon density will likely stay the same 
(unless there is diffusion of stars from the center outward) . 

If all dwarf galaxies are structurally similar (which means that the integral in eq. |||] is the 
same for all dwarfs), and in the core almost all the gas is transformed into stars {pb,c ~ P*,c, this 
assumption is addressed below), then equation (^ reduces to the power-law dependence between 
the total luminosity of a galaxy and its central luminosity density Ic, 

LJtI^I^ (3) 

at t = tf. However, this equation also holds at the present time if the star formation at later times 
{t > tf) is negligible or can be corrected for. 

It is important to point out here, that it is the luminosity density (i.e. the stellar mass 
density) and not the total mass density that matters. The dark matter density in the Local Group 
dwarfs varies enormously and is not correlated to the stellar density, as can be easily seen from the 
total mass-to- light ratios (Table 4 of Mateo 1998). Since the galaxies included in Table |^ formed 
most of their stars more than 10 Gyr ago, it can be presumed that the observed core luminosity 
densities of the dwarfs are representative of the core luminosity densities at t/ ~ 3 Gyr. This 
cannot be said about the dark matter, however. In the last 5 Gyr or so the Local Group dwarfs 
inhabited a region with the mean density about hundred times higher than the mean density 
of the universe, and therefore they had enough opportunities to accrete large quantities of dark 
matter. Thus, their dark matter contents is representative of their whole evolutionary history, and 
not of only the first 3 Gyr of their life. 

The data from Table || can be used to investigate the relationship (^) in the Local Group 
dwarfs. It is natural to take tf as the time 10 Gyr ago, corresponding to the sharp drop in the star 
formation rate in all galaxies in Table ||. For a 13 Gyr old universe this would imply tf = 3 Gyr. 
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Fig. 1 . — The corrected V-band luminosity per cube of the core radius versus the central luminosity density 
for the Local Group dwarfs from Table |l|. The dotted line shows the power-law fit with the slope fixed to 
3/2. 

However, even if the galaxies in Table |l] formed most of their stars more than 10 Gyr ago, 
they still had additional star formation at later times. These "recent" stars need to be excluded 
from the analysis. It is straightforward to estimate the fraction />io of all stars formed more than 
10 Gyr ago, from Fig. 8 of Mateo (1998) and Fig. 5-6 of Grebel (1998). The quantity />io is listed 
in the last column of Table ||. Thus, instead of L* and Ic from equation |3|, 1 use = L^f^iQ and 
Ic = Icf >io instead. Figure |l| shows the relationship (^) for the data from Table 1 together with a 
power-law fit with the slope n = 3/2. 

However, Fig. |l] is a bad way to represent the data because the errors along the two axes are 
not independent - both quantities include the large errors due to uncertainties in the distances 
and the core radii. It is better to represent the relationship (^) as 



y = 



-'c'c 



n-1 



(4) 
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In this form the left hand side is dimensionless, and thus independent of the distance, whereas the 
errors on Ic are dominated by the error in the distance, because the central luminosity density 
is inversely proportional to the distance cubed. In this form the errors on y and are not 
strongly correlated, and I can measure the statistical significance of the correlation between the 
two quantities in equation (Q). 

Given the the errors on the observed quantities: the relative photometric error 5V in the total 
luminosity, the relative error in the measured angular core radius Stc-, the relative distance error 
5D, and the relative photometric error in the observed central surface brightness 5Tic, the errors 
on y and 1^ are calculated as = SV"^ + A5rl + 5Til, and 5ll = 5Til + 5rl + 95 D'^. 

However, my correction to the total luminosity and the central luminosity density is, 
admittedly, uncertain, because the star formation histories are not known too accurately, and thus 
the quantity />io is somewhat uncertain as well. To account for this uncertainty, I increase the 
errors on Ic a factor of l//>io, in which case the contribution of dwarfs with a large correction to 
the total is reduced by a factor of />io- 

The best power-law fit to the data listed in Table |^ is given by the following formula: 

In this form the errors on the amplitude and the slope are uncorrelated. If, instead, I use 
uncorrected errors, the reduction of the fit error is minor: the power-law index is 0.52 it 0.10, and 
the logarithm of the amplitude is 1.20 it 0.06, which demonstrates that my correction procedure 
is not too important. I can also fit the original data, without correcting the total luminosity and 
the central luminosity density by a factor />io. In the latter case I obtain an equally good fit: the 
power-law index is 0.51 ± 0.10, and the logarithm of the amplitude is 1.16 it 0.06. The reason why 
the original data work equally well is that the fit is dominated by the data points with the smallest 
error-bars, and those also have the smallest correction. Whether this is a mere coincidence, or has 
a physical meaning, can only be speculated at the moment. 

All three fits are highly significant statistically, which indicates that the error bars are not 
underestimated . 

The result (|5|) is non-trivial. It shows that the average luminosity density of dwarf galaxies 
was proportional to the 3/2 power of their central luminosity density in the first 3Gyr. The 
central luminosity density, as I argue below, is representative of their original gas density, since the 
star formation at the center is efficient enough to be able to convert all their gas into stars. The 
average luminosity density (i.e. the total luminosity) of the dwarfs is however dominated by the 
outer regions (as can be seen from Table |l| since y ^ 1). In the outer regions the star formation is 
inefficient, and thus the total luminosity is proportional to the gas density (and thus the central 
luminosity density) to the 3/2 power (or about that), and not to the central luminosity density to 
the first power, which would be the case if the star formation was efficient everywhere and all the 
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gas in the dwarfs was converted into stars. 

It is remarkable that the observed slope is so close to the orthodox Schmidt law power-law 
index of 3/2. Therefore, it is worthwhile to investigate this special case in more detail. The 
star formation rate in this case can be rewritten in the following form, fixing the coefficient of 
proportionality in the Schmidt law: 

aZ Cdyn 

where a dimensionless quantity e measures the efficiency of star formation, and the dynamical 
time idyn is defined as in Binney & Tremaine (1987), tdyn = {3TT/[l6Gpg]y/^ . Note, that I define 
the dynamical time as depending on the gas density, and not on the total density. If one desires to 
define the dynamical time in terms of the total density, e simply needs to be decreased by a factor 

of {no/nb)^/\ 

Integrating equation (|6|) over time and space, and again invoking At^s to hide the consumption 
or infall of gas, I obtain the following relationship between the total luminosity of a galaxy and 
the core baryon density: 

= 3.37 X 10-^ cm3/2 g-i/2 s-^ttT-^ At.sepl^Jr^Cy^, (7) 

where T is the mass-to-light ratio, and C„ = {Pb/ Pb,c)^ f"^ dr /r^. For the King model with the 
concentration parameter 2.5, C3/2 = 0.6. Since, as I argue below, the central luminosity density in 
the dwarfs is representative of the central baryon density (i.e. almost all of the gas is consumed 
at the center), the King profile is a good approximation to the total baryon density profile near 
the center. The Navarro, Prenk, & White (1997) density profile, motivated by cosmological 
simulations, also gives C3/2 = 0.6, which indicates that this number is not very sensitive to the 
details of the density profile. 

Finally, I can relate the core baryon density to the core luminosity density by ph^c = '^Ic/ f*, 
where Z^, is the fraction of baryons turned into stars at the center. 

Comparing equation (|5|) taken with the power-law index 3/2 with equation (|7|), I obtain the 
expression for the value of the star formation efficiency e: 

e = 0.040 X 10^° °^At^iT4 (8) 
where At2 = Ateff/2Gyr, and T4 = T/{AMq/Lq). 

Using equations ^ and (|8|), I can obtain the following expression for the star formation time 
T^,^c = p*^c/{dp*^c/dt) at the center of a dwarf galaxy: 

r,,, = 1.5 X 10^ GyrAt2/2"'/'(l - /*)"'/', 

where I2 = Ic/{0.01Lq/ pc^). Thus, for all dwarfs with /c ^ 2 x IO^^Lq/pc'^ (which are all 
galaxies listed in Table |l| with the exception of the Sextans and Ursa Minor dwarf spheroidals) 
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■^*,c ^ tf until most of the gas is consumed and becomes close to unity (an average la error 
on Ic corresponds to about 40% error on so = 0.6 and /* = 1 cannot be distinguished on 
the basis of the current data). Thus, assuming = 1 (which is certainly true for the most dense 
dwarfs with Ic"^ 2 x lO~^L0/pc^), I obtain the final expression for the average efficiency of star 
formation in the first 3 Gyr: 

eAt2 = 0.040 X 10±°°''T~^/^. (9) 

While T4 is unlikely to be different from unity, the last remaining factor At2 is less certain. 
It cannot be larger than 1.5, since the gas density cannot exceed the total baryon density. In 

— 1/2 

particular, e > O.O22T4 at 3 o". However, At2 can be much smaller than 1 (giving e ^ 0.04) 
if most of the gas in the dwarfs was consumed (and never replenished) in a short time interval 
well before the universe turned 3 Gyr (it is not clear whether the existing data rule out such a 
scenario). 



3. Conclusions 

I use a sample of Local Group dwarf galaxies selected by the criterion that they form most of 
their stars more than about 10 Gyr ago to measure the star formation law in the first 3 Gyr of the 
life of the universe (assuming the age of the universe is 13 Gyr). I find that the star formation rate 
on scales comparable to the sizes of dwarf galaxies (100 — 600 pc) is a power-law function of the gas 
density, with the power-law index 1.54 it 0.11. If the index is fixed at n = 3/2, the star formation 
rate is given by equations (|6|) with an efficiency of about 4% if the gas is consumed continuously. 
The efficiency cannot be smaller than 2.2% but can be higher if the gas was consumed well before 
the universe turned 3 Gyr. 

The efficiency is not a strong function of scale for scales in the range from about 100 to 600 
parsecs (if it were, the straight power-law would not be a good fit in Fig. 

Admittedly, the time of the drop in the star formation rate in the Local Group dwarfs is 
not very well determined, and is, perhaps, uncertain to about one or two billion years. Thus, 
the "simultaneous" drop in star formation rates in the majority of the dwarfs is perhaps not 
sufficiently simultaneous to preclude the possibility that it was supernova driven winds rather 
than reheating of the universe (or local reionization) , which was responsible for the observed drop. 
However, in this case it would be difficult to explain the observed correlation, since the central 
regions of the dwarfs galaxies can be expected to be affected most by the supernova activity, thus 
spoiling the beautiful correlation shown in Fig. |l|. Therefore, the present result indicates that the 
supernova driven winds play only a modest role in the early evolution of low-mass galaxies. 

Curiously, equation ^ can be used, in principle, as a novel distance indicator, since only the 
X-axis depends on the distance. However, it is likely to not be competitive with other indicators 
given the amount of information needed to measure Ly and Ic- 
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The result presented here can be incorrect if the structural form of the Local Group dwarfs 
changes monotonically with their luminosity as to mimic the observed relationship (^). I would 
like to emphasize that, while this possibility cannot be completely excluded, there is currently 
no observational evidence that supports it, and this assumption also contradicts the modern 
theoretical view on the structure of galaxies. 



I am grateful to John Bally, Ed Bertschinger, Brad Gibson, Oleg Gnedin, Phil Maloney, Mike 
Shull, John Stocke, and the anonymous referee for useful comments and discussions. 



REFERENCES 

Binney, J., & Tremaine, S. 1987, Galactic Dynamics (Princeton: Princeton University Press), 37 
Contardo, G., Steinmetz, M., & Fritze-von Alvensleben, U. 1998, ApJ, 507, 497 
Gnedin, N. Y. 2000, ApJ, submitted ( |astro-ph/0002151| ) 
Grebel, E. 1998, lAUS, 192, 1 

Hernandez, X., Gilmore, G., & Valls-Gabaud, D. 2000, MNRAS, in press ( ^stro-ph/0001337| ) 

Mateo, M. 1998, ARA&A, 36, 435 

McWilliam, A., & Rich, R. M. 1994, ApJS, 91, 749 

Navarro, J., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 

Navarro, J., & Steinmetz, M. 1997, ApJ, 478, 13 

Quinn, T., Katz, N., & Efstathiou, G. 1996, ApJ, 278, 49P 

Ricotti, M., Gnedin, N. Y., & Shuh, J. M. 2000, ApJ, in press ( |astro-ph/9906413| ) 

Schaye, J., Theuns, T., Ranch, M., Efstathiou, G., & Sargent, W. L. W. 2000, MNRAS, submitted 
(|astro-ph/9912432|) 

Thoul, A. A., & Weinberg, D. H. 1996, ApJ, 465, 608 

van den Bergh, S. 1994, ApJ, 428, 617 

van den Bergh, S. 1999, aRev., 9, 273 

Weinberg, D. H., Hernquist, L., & Katz, N. 1997, ApJ, 477, 8 



This preprint was prepared with the AAS lATp^X macros v4.0. 



